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calculations are in good agreement with experimental data. 
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1 Introduction 



The reaction 7 Be(p,7) 8 B plays an important role in the so-called ppIII chain 
in the hydrogen burning of main-sequence stars. A knowledge of the reaction 
rate is essential to determine the branching ratios between the ppl, ppll and 
ppIII chains. The magnitude of the reaction cross section is of special interest 
for the solar neutrino problem. The reason for this is that in the 37 CI neutrino 
experiment 77%, in the Kamiokande II experiment 100% and in the 
gallium experiments @, |J 11% of the detected neutrino flux originate from 
the high-energy neutrinos emitted in the ppIII chain ||. Therefore, the 
reaction rate of the above process determines the high-energy solar neutrino 
flux. 

The reaction considered in this paper has been measured by various au- 
thors at subCoulomb energies. The most recent data have been obtained in 
the energy range E c _ m , = 117 — 1230 keV 0. 

The dominance of the direct interaction (DI) mechanism and the validity 
of the description with the potential-model approach below the Coulomb 
barrier has been established in many light-ion reactions (0 and references 
therein). In this work we apply this approach to the above reaction. In the 
next section we introduce the direct capture model. In Section 3 the result 
for the astrophysical ^-factor is given and compared to experimental data. 
Finally, in Section 4 the results are summarized. 
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2 Potential Model Approach 

Potential models are based on the description of the dynamics of the reaction 
by a Schrodinger equation with local optical potentials in the entrance and/or 
exit channels. Such models are the "Distorted Wave Born Approximation" 
(DWBA) 1, |, g for transfer or the "Direct Capture" model (DC) 0, 0, 



for capture reactions. 
In first order perturbation theory the expression for the differential cross 
section of a direct capture reaction is ]14 : 



dO, ~ {He) U c) UJ 2I A + 1 2S a + 1 m ± Mb(t 1 1maM « M »>° 1 ' 

(1) 

Here I a {Ma) and S a (M a ) are the spins (their projections on the z-axis) 
of target and projectile, and a is the polarization of the electromagnetic 
radiation (cr = ±1). The wave numbers of the emitted 7-rays and of the 
asymptotic relative wave function in the entrance channel are denoted by A; 7 
and k a , and m is the reduced mass. We couple the angular momenta in the 
spin-orbit representation: 

la + S a =J a , (2) 

h + S a = j b , (3) 
Ia + Ib =Jb ■ (4) 



The coupling and notation are the same as in |L4 



The transition amplitudes TM A M a M B ,o are expanded in terms of rotation 
matrix elements d} a (8) with the electromagnetic multipole A (A = El, E2, 
Ml, ...), 

TM A M a M B ,(T = ^2T MAMaMB;(T d Sa (9) , (5) 

A 

where 5 = M A + M a — Mb and 9 is the angle between k a and k 7 . 
The transition amplitudes are proportional to the radial integrals 

T " « /,:/.„., =Jdr U hjblB (r) O x (r) Xlaja (r) , (6) 

where Ui b j b i B (r) and Xi a j a ( r ) are the radial parts of the bound state wave func- 
tion and the distorted wave function in the entrance channel, respectively. 
The functions O x (r) are the radial parts of the electromagnetic multipole 
operators, which were taken in their approximated form, 

C M1 (r)~l, C E1 (r)~r, C E2 (r) ~ r 2 . (7) 

Folding potentials 

We solve the radial integral @ numerically using single folding potentials |15 
for the bound state potential and the optical potential, 

V(r) = A J dr APA (r A ) t(E, p A , |s = r - r A |) , (8) 

which were assumed to be real. Here r is the separation of the centers of 
mass of the bound or colliding particles, p A is the nucleon density of the 
target and A is a normalization constant. For the effective nucleon-nucleon 



interaction t in the entrance and exit channel we chose the density dependent 



form of the M3Y interaction 16 



The nucleon density pa is derived from the charge density distribution 
p p (r) of the target A with the assumption, that the distribution in the nucleus 
is the same for protons and neutrons: p n = (N/Z)p p [§, p. 474]. This leads 
to 

PA = P P + (N/Z) Pp = (A/Z) p p . (9) 
The density Pa(j') satisfies 

4tt J p A (r)r 2 dr = A, (10) 

The normalization constant A of the folding potential (A sc for the potential 
in the entrance channel and Ab for the bound state potential) accounts for the 
interplay of the Pauli principle and distortion and breakup effects. One of the 
advantages of the folding procedure lies in the fact that no open geometrical 
parameters exist, since Ab is adjusted to reproduce the experimentally known 
binding energies and A sc can be adjusted to reproduce scattering phase shifts 
(if they are available) or the energies of resonant states (which is done for 
7 Be(p,7) 8 B). Therefore the form of the optical and bound state potentials is 
determined uniquely. 
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3 Results 

7 Be(p, 7 ) 8 B 

The 7 Be(p,7) 8 B reaction is of crucial importance since it leads to the high- 
energy 8 B neutrinos. Unfortunately, there are still significant experimental 
uncertainties in the low-energy cross section for this reaction. The reaction 
was analyzed before in the DC model by using phenomenological poten- 
tials |17], jl8| . Furthermore, this reaction was also analyzed with the res- 
onating group method (RGM) |L9] and the generator coordinate method 
(GCM) [fJU], ^TJ. The numerical calculations in the DI model in this work 



were performed using the direct capture code TEDCA [52]. 
The folding potential for 7 Be+p 

Because of the instability of 7 Be there are no experimental charge distribu- 
tions available to compute the nucleon density with Eq. @ p3| . Therefore 
we chose the same distribution as that of 7 Li, which should give a suitable 
approximation. Although the experimental charge distributions of 7 Li and 
7 Be are somewhat different, the resulting potential should be quite similar. 
We conclude this from calculations comparing the folding potentials of 11 B+p 
and 12 C+p. 

The folding potential for 7 Be+p was calculated using a harmonic oscillator 
model for the 7 Li charge distribution with parameters a = 1.77 fm and 
a = 0.327 [|||. The depth of the pure folding potential is 43 MeV (A = 1), 



6 



the volume integral per nucleoli is 660 MeV fm 3 . 
Calculation of 'Be(p,7) 8 B without spin-orbit coupling 

All contributions to the transition matrix lead to the ground state of 8 B 
with spin J w = 2 + . We use the spectroscopic factors Si P3/2 = 0.977 and 
<Si Pl/2 = 0.0561 from the shell model calculations of |2H}| . The DC cross 



section is a sum over the two final state configurations (jj, = 1/2, 3/2), 

ct dc = E^Xi^J! > ( n ) 

hjb 

where C 2 is the isospin Clebsch-Gordan coefficient. 

A strength factor of Ab = 1.041 yields the correct separation energy of 
138 keV of the bound proton. The bound state wave function is very diluted 
due to the very low binding energy. Therefore we had to calculate radial 
contributions to the transition matrix up to 200 fm. 

In Fig. 1 the astrophysical 5-factor is shown. The A-parameter of the 
scattered wave A sc = 0.969 was determined with the condition that the res- 
onance in the p-wave has to be at an energy of 632 keV, which leads to 
a resonant Ml transition. This value of the normalization parameter de- 
termines the potential uniquely and was used for all partial waves in the 
entrance channel. 

In Fig. 2 the radial integrand (6) of the s-wave is plotted at two differ- 
ent energies (E c . m , = 15 keV and i? c . m . = 1 MeV). At low energies there 
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are contributions to the integral from very far outside the nucleus. Those 
contributions can be neglected only for distances larger than 200 fm. 

As one can see in Fig. 1, the calculated S'-factor of 7 Be(p,7) 8 B agrees 



with the experiments || |25|, £6| quite nicely in absolute value over a wide 
range of energies (100 keV — 2 MeV), except in the resonance energy region 
near -E c .m. = 632 keV. In particular there is also a very good agreement 
with the mean experimental value of 5*(0) that is used by Bahcall || in the 
standard solar model. A comparison of the parametrization of the potential 
model calculation with these experimental values is given in Tab. 1. The 
data of [HI in Fig. 1 were renormalized as suggested by |27|. 

The experimental width of the l + -resonance at 632 keV is not reproduced 
so well. The calculated width of 100 keV is too large by a factor of 2.5. A 
possible explanation for this is the neglect of configuration mixing and other 
channels, such as the channel p + 7 Ben/ 2 )-, which could only be taken into 
account in a coupled channel calculation. 

Calculation with spin-orbit coupling in the entrance channel 

Currently another resonance is discussed, which should correspond to a 
J w = 1 + state of 8 B with an excitation energy E x pa 1.5 MeV. This state 
is derived from comparison with the excitation spectra of the mirror nucleus 
8 Li |H (If : £ x = 0.9808 MeV, l£ : E x = 2.255 MeV). Up to now, no 
experimental evidence for the 1^ state in 8 B is available. The lj state is 
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predicted by the GCM |20| and also by the potential model, if one takes 



spin-orbit coupling into account. The two states are assumed to have the 
configurations 7 Be+p(lp3/2) (1^) and 7 Be+p(lpi/2) (1^)- 

We are able to include this doublet in our calculation by employing an 
additional spin-orbit potential in the entrance channel. For this spin-orbit 



potential we used the usual parametrization |£9| , with the Woods-Saxon form 
factor replaced by the form factor of the folding potential. The strength 
A sc = 0.94, together with the spin-orbit term (V so = 2.44 Me V), reproduces 
the correct energy separation of the two l + -resonances at -E c . m . = 632 keV 
and -Ec.m. = 1.4 MeV. Fig. 3 displays a plot of the contributions of different 
partial waves and multipolarities. The influence of the J n = l^-resonance 
on the S'-factor is negligible, especially for thermonuclear energies. 

Calculation with spin-orbit coupling in both channels 

According to the shell model calculations of |2"4]|, the 1^ level in 8 B has a 



better p 3 / 2 structure (<Si P3/2 = 0.3215 and <5>i Pl/2 = 0.1240). Nevertheless, 
we also want to consider the description of the ground and first excited 
(unbound) state of 8 B as a p-wave spin-doublet. The two states should then 
have the configurations 7 Be+p(lp3/2) (2 + ) and 7 Be+p(lpi/ 2 ) (1^). 

The form factor of the spin-orbit potential stays the same as for the 
calculation with spin-orbit coupling in the entrance channel. The spin-orbit 
strength V^ Q = 2.01 MeV was adjusted to describe the correct energy splitting 
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of the two corresponding states in 8 Li. The so determined spin-orbit potential 
is able to reproduce with Ab = A sc = 1.01733 both the binding energy of the 
ground state and the resonance in the entrance channel at 0.632 MeV. The 
potentials for the three assumptions (no spin-orbit potential, spin-orbit in 
the entrance channel and spin-orbit in both channels) are shown in Tab. 2. 

Fig. 4 displays the ^-factor calculated with the additional spin-orbit po- 
tential in both channels. Now the height of the resonance at 632 keV is 
reduced, compared to the previous calculations. At this energy, the dom- 
inating Ml-contribution is formed by the overlap of the quasi-bound and 
bound state (Eq. (6)), which are now almost orthogonal. This is due to 
the assumption of a spin-dublet with a spin-orbit potential. In Fig. 5, the 
radial wave function of the bound state (full line) is shown together with 
the scattering wave functions of the p wave at resonance energy (632 keV) 
for the case of zero spin-orbit potential (dashed line) and spin-orbit potential 
(dashed dotted line). The overlap of scattering wave and bound state reaches 
its maximum at 3 fm in both cases, but is much smaller in its absolute value 
for the case of spin-orbit potential. 

As a consequence, the description of the experimental S'-factor by the 
potential model gets worse in the resonance energy region (Fig (4)). The 
assumption of a spin-dublet does not lead to a good description of the data. 
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7 Li(n, 7 ) 8 Li 

We now want to calculate the mirror reaction 7 Li(n,7) 8 Li with the assumption 
of spin-orbit coupling in the entrance channel only, because this seems to be 
the most realistic potential. Since the determination of the folding potential 
is not depending on isospin, the 7 Li+n-potential is the same as the 7 Be+p- 
potential of the last section. Therfore the optical potential for the entrance 
channel can be taken from the analog calculation of 7 Be(p,7) 8 B (see also 
Tab. 2). Also the shell model spectroscopic factors are the same for 8 Li and 
8 B. 

We consider the capture into the ground state (J 71 = 2 + ) and the first 
excited state (J 71 = 1+ E x = 0.9808 MeV) of 8 Li. The correct binding 
energies of these states are reproduced with the parameters Ab = 1.0361 for 
the ground state and Ab = 0.9627 for the first excited state of 8 Li. 

The result of the calculation is shown in Fig. 6. The experimental data 
were taken from [|30| (filled circles) and [[H]]. The open circles and triangles 
result from two different normalization procedures. We did not attempt to 
reproduce the resonance at E c , m . = 0.222 MeV, which corresponds to a 3 + 



state in Li, since this state has a bad Li+n structure 32 
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4 Summary 

We have shown that the reactions 7 Be(p,7) 8 B and 7 Li(n,7) 8 Li can be de- 
scribed by a potential model at thermonuclear energies. The theoretical 
values of the low-energy astrophysical .S-factor for the different assumptions 
concerning the single-particle configurations of 8 B are in good agreement 
with the mean experimental values of ||. It can be concluded that the di- 
rect interaction mechanism is dominant at energies well below the Coulomb 
barrier. 

The resonance width and strength of the experimentally known resonance 
at i? c . m . = 632 keV are both overestimated by the potential model, if one 
uses no spin-orbit potential or a spin-orbit potential in the entrance channel 
only. In these cases no special assumptions concerning the single-particle 
configurations of the final bound state in 8 B are made. If one assumes the 
configurations 7 Be+p(lp3/2) (2 + ) and 7 Be+p(lpi/2) (1^) for this state and 
introduces an additional spin-orbit potential to describe the bound and first 
unbound state simultaneously, the DC cross section in the energy region 
of the resonance is strongly reduced (Fig. 4) and underestimates the data, 
due to the overlap of relatively orthogonal wavefunctions. But in the non- 
resonant energy region the reproduction of the cross section data is excellent 
in all cases. 

The reaction 7 Li(n,7) 8 Li was calculated with the potential in the entrance 
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channel taken from 7 Be(p,7) 8 B, since this seems to be the most realistic and 
therefore favorable potential ansatz. The agreement with the reaction data 
is good, as can be seen in Fig. 6. 

The progress of applying the potential model combined with the folding 
procedure lies mainly in the fact that no parameter has to be adjusted to the 
reaction data. This is neither the case in the phenomenological fits of Barker 



et al. |27J nor Kim et al. [14]], where potential parameters have to be fitted 
to the experimental reaction data. 

In our model the strengths of the folding potentials are adjusted to re- 
produce the energies of the bound and quasi-bound states. This procedure 



is similar to microscopic models P0|| , where the correct binding energies are 
ensured by adjusting a parameter in the nucleon-nucleon interaction. 
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Figure captions: 

Fig. 1: Astrophysical S'-factor for the reaction 7 Be(p,7) 8 B. Full curve: poten- 
tial model calculation (no spin-orbit coupling). The data points are from || 
(triangles), |25[ (circles) and [^(J (squares). 

Fig. 2: 7 Be(p,7) 8 B: Radial integrand of Eq. (6) for two different energies in 
the entrance channel (s-wave only). 

Fig. 3: Contributions of different partial waves and multipolarities to the 
cross section of 7 Be(p,7) 8 B. Here an additional spin-orbit potential in the 
entrance channel was used which influences the p-wave contribution near 
E c . m . = 1.4 MeV. 

Fig. 4: Astrophysical S'-factor of 7 Be(p,7) 8 B. Full curve: potential model 
calculation (spin-orbit coupling in both channels). 

Fig. 5: Radial wave functions of the 7 Be+p bound state (full line) and scat- 
tering wave functions of the p wave at resonance energy (0.632 MeV) for zero 
spin-orbit potential (dashed line) and spin-orbit potential in both channels 
(dashed dotted line). 

Fig. 6: Total cross section for the reaction 7 Li(n,7) 8 Li. The experimental 
data are from f3(J (filled circles) and f3TJ (open circles and triangles). 



Table captions: 

Tab. 1: Comparison of the results of the potential model calculation of 
7 Be(p,7) 8 B for the three different assumptions (no spin-orbit coupling, spin- 
orbit in the entrance channel and spin-orbit in both channels) with the ex- 
perimental values of the astrophysical ^-factor and its derivative as given 
in §. 

Tab. 2: Comparison of the potentials used in the calculation of 7 Be(p,7) 8 B 
and 7 Li(n,7) 8 Li. 
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